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ABSTRACT 

We solve for the velocity fields of momentum-conserving supershells driven from galaxy 
centres by steady winds from supermassive black holes or nuclear star clusters (central 
massive objects: CMOs). We look for the critical CMO mass that allows such a shell 
to escape from its host galaxy. In the case that the host galaxy dark matter halo is 
a singular isothermal sphere, we find that the critical CMO mass derived by King, 
which scales with the halo velocity dispersion as Merit oc cr'*, is necessary, but not by 
itself sufficient, to drive shells to large radii in the halo. Furthermore, a CMO mass 
at least 3 times the King value is required to drive the shell to the escape speed of 
the halo. In the case of CMOs embedded in protogalaxies with non-isothermal dark 
matter haloes, which we treat here for the first time, we find a critical CMO mass 
that is sufficient to drive any shell (under a steady wind) to escape any galaxy with a 
peaked circular speed profile. In the limit of large halo mass, relevant to real galaxies, 
this critical CMO mass depends only on the value of the peak circular speed of the 
halo, scaling as Merit V^^^^. Our results therefore relate to observational scalings 
between black hole mass and asymptotic circular speed in galaxy spheroids. They 
also suggest a natural way of extending analyses of M-a relations for black holes in 
massive bulges, to include similar relations for nuclear clusters in lower-mass and disc 
galaxies. 
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1 INTRODUCTION 

Most early- type galaxies and bulges with M > 10^° Mq har- 
bour a supermassive black hole (SMBH) at their centre (Ko- 
rmendy & Richstone 1995); while observations with the Hub- 
ble Space Telescope have revealed the presence of massive 
nuclear star clusters (NCs) in the majority of less massive 
galaxies (both early- and late-type: Phillips et al. 1996; Car- 
oUo et al. 1997; Boker et al. 2002; Cote et al. 2006, 2007). 
The properties of these central massive objects (CMOs) cor- 
relate tightly with properties of their host galaxies, perhaps 
most notably in terms of CMO mass, M;,^^ , versus (bulge) 
stellar velocity dispersion, a: oc , with a; ~ 4 (for 

SMBHs, see, e.g., Ferrarese & Merritt 2000, Gebhardt et 
al. 2000, Tremaine et al. 2002, Ferrarese & Ford 2005, or 
Giiltekin et al. 2009; for NCs, see Ferrarese et al. 2006; also 
relevant are Wehner & Harris 2006 and Rossa et al. 2006). 
Though essentially parallel, there is an offset between the 
M-cr relations of NCs and SMBHs, in the sense that the 
NC masses in intermediate- and low-mass galaxies tend to 
be ~ lOx larger than if they followed a simple extrapola- 
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tion of the SMBH Ad-a relation for higher-mass spheroids 
(Ferrarese et al. 2006; see also McLaughlin et al. 2006). 

Recently, Volonteri, Natarajan, & Giiltekin (2011; 
cf. Ferrarese 2002) have argued that galaxies/bulges contain- 
ing SMBHs also show a correlation, of the form Mbh oc 
with y ~ 4, between black hole mass and the "asymptotic" 
circular speed 14 at large radii where dark matter is ex- 
pected to dominate the total galaxy mass. There is some 
debate (e.g., see Ho 2007; Kormendy & Bender 2011) over 
how the stellar velocity dispersions in the M-a relation, 
which are measured inside a fraction of the bulge effective 
radius, connect empirically to asymptotic circular speeds, 
which normally refer to many times the effective radius de- 
fined by stars. This can be a difficult question (with a model- 
dependent answer), especially in "hot" stellar systems where 
circular speeds — that is, Vc{r) = GM{r)/r — are not ob- 
served simply as net rotation. However, the existence of 
some kind of connection, and at least the possibility of an 
Mbh-Vc relation in addition to Mbh-o", is clear in principle: 
The stellar velocity dispersion at any radius in a dark-matter 
dominated galaxy depends on the dark matter distribution, 
which is precisely what Vc probes at large radii. 

Self-regulated feedback from growing CMOs is thought 
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to play a key role in establishing the M-a relation and as- 
sociated scalings. Though through different mechanisms, ei- 
ther an NC or an SMBH will drive an outflow, which sweeps 
the ambient gas in a protogalaxy into a shell that, at least 
initially, is able to cool rapidly and is therefore momentum- 
driven (King 2003; McLaughlin et al. 2006; see also ig] be- 
low). There is then a critical CMO mass above which the 
outwards force of the wind on the shell may overcome the 
inwards gravitational pull of the CMO plus the dark matter 
halo of the parent galaxy. 

The only case that has been considered in detail analyt- 
ically is that of a steady wind, in which the CMO mass (and 
associated wind thrust) is constant throughout the motion 
of the shell (Silk & Rees 1998; Fabian 1999; King 2003, 2005, 
2010a; Murray et al. 2005; McLaughlin et al. 2006, Silk & 
Nusser 2010). In this case, and assuming a halo modelled as 
a singular isothermal sphere (SIS) with velocity dispersion 
CTo, King (2005) found a critical CMO mass of 
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(see also Mclaughlin et al. 2006; Murray et al. 2005). In 
this expression, k — 0.398 cm^g~^ is the electron scatter- 
ing opacity; fo is an average gas mass fraction (« 0.2, so 
fo.2 = /0/O.2); and (T200 ~ cro/200 km s~^. The parameter 
A is related to the feedback efficiency for each type of CMO; 
it has a value X ^ 1 for SMBHs, and A ^ 0.05 for NCs 
(McLaughlin et al. 2006). Once a CMO in an isothermal 
halo with a given cto has grown to at least the mass in equa- 
tion IT}, the CMO wind may drive a momentum-conserving 
shell with coasting speed u > at arbitrarily large radii in 
the galaxy. This then admits the possibility of a blow-out 
clearing the galaxy of any remaining ambient gas, choking 
off further star formation and CMO growth, and locking in 
an Mp„Q-cr relation. 

As a momentum-driven shell moves outwards from a 
CMO, gas cooling times increase and a switch to an energy- 
driven phase is expected, at which point the shell can accel- 
erate to escape the galaxy (King 2003). Momentum-driving 
may then need only push a shell out to where the switch 
to energy- driving occurs; and this can be done with a CMO 
less massive than the Merit in equation IT}, which is neces- 
sary for momentum-driving to arbitrarily large radii. This 
suggests that equation IT} may actually predict an upper 
limit for observed M-a relations; and indeed, the equation 
lies above current best fits to data by factors of a few. 

Distributed star formation in a protogalaxy bulge is ex- 
pected to provide additional momentum input to the feed- 
back (Murray et al. 2005; Power et al. 2011). This would also 
reduce the CMO mass required for the feedback to escape, 
again suggesting that the M^mq-ct relation in equation |T]) 
is an upper limit. 

Silk & Nusser (2010) have shown that, in a truncated 
isothermal sphere specifically, a momentum-conserving shell 
driven solely by a steady black-hole wind can reach large 
radius with fast enough speed to escape directly (that is, 
with V > 2(To), only if the SMBH mass is at least a few 
times larger than the critical value in equation |l} (which is 
necessary just to have d > at large r). This would put the 
predicted normalization of an M^^jq-ct relation above the 
observed normalization by a full order of magnitude. Silk & 
Nusser argue from this that the real key to a feedback origin 



formation that is triggered by the outfiow from a CMO. 
However, Power et al. (2011) counter that a switch from 
momentum- to energy-driving of the CMO feedback is still 
inevitable and will alleviate some of the difficulty identified 
by Silk & Nusser. 

In this paper, we investigate how this basic feedback sce- 
nario for M-a relations depends on the simplifying assump- 
tion that dark matter haloes are SISs. We analyze aspects of 
the dynamics of supershells in spherical but non-isothermal 
haloes, while retaining some other simplifying assumptions 
(steady winds and purely momentum-driven shells) in com- 
mon with previous work. 

Our main result is a generalization of the critical CMO 
mass that suffices to blow momentum-driven feedback en- 
tirely out of any realistic, non-isothermal dark matter halo 
that has a well-defined maximum in its circular speed pro- 
file, V^{r) — GM^j^ (x)/'''- For large halo masses, this critical 
CMO mass tends to the limiting value. 
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where 14, pk is the peak value of the circular speed. 

In a SIS, which has a constant Vc = \/2cro, our new 
equation |5)) clearly reduces to equation |T}. However — as 
we discuss in detail in !j3]and i ]4. II below — in an isothermal 
halo this Merit is necessary hut not sufficient, in general, to 
guarantee the escape of a momentum-driven CMO wind. By 
contrast, in the more realistic, non-isothermal cases that we 
consider, equation |2]) gives the M^^^ that is sufficient for 
the escape of any such feedback. 

Any momentum-conserving shell driven by a steady 
wind from a CMO with the mass in equation |2| will even- 
tually accelerate at large radii and exceed the escape speed 
of any non- isothermal halo with a peaked Vc{r) profile, even 
without a possible change to energy-driving, additional mo- 
mentum feedback from star formation or growth of the CMO 
(none of which we include in our analysis). Thus, the objec- 
tion of Silk & Nusser (2010) to equation |l} as the basis for 
observed M-a relations applies only if dark matter haloes 
are strictly isothermal. 

Equation |2} defines the "characteristic" velocity dis- 
persion that needs to be considered when interpreting ob- 
served M(-,^jo-o" relations in non-isothermal galaxies: ao = 
Ve,pk/%/2. It also gives the first direct, quantitative pre- 
diction of an M^j^Q-Vc relation such as that discussed by 
Volonteri et al. (2011). The result may still be an upper limit 
to observed relations since we do not consider any transition 
to energy-conserving feedback, nor any sources of feedback 
other than steady CMO winds, in this paper. 

We begin in by looking at the general equation of 
motion of a momentum-driven shell as it moves out into a 
gaseous protogalaxy. In !j3]we develop, in more detail than 
before, the case of the SIS. In H4.ll we analyze the motion of 
a momentum-driven shell in a general, non-isothermal halo 
with a peaked circular-speed curve, and derive equation |2]). 
In the rest of 331 we illustrate our general results using three 
particular dark-matter halo models as examples (those of 
Hernquist 1990; Navarro, Frenk, & White 1996, 1997; and 
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Dehnen & McLaughlin 2005). In Sj5]we summarize the paper 
and give a brief discussion. 



2 EQUATION OF MOTION 

An SMBH accreting at near- or super-Eddington rates in a 
gaseous protogalaxy is expected to drive a fast wind back 
into the galaxy (King & Pounds 2003), with quasi-spherical 
(i.e., not highly coUimated) geometries indicated by obser- 
vations of strong outflows from local AGN (e.g., Tombesi 
et al. 2010). Similarly, the combined winds and supernovae 
from massive stars in a very young (still forming) NC will 
drive a superwind into its host protogalaxy. In a spherical 
approximation to either case, the wind sweeps up the sur- 
rounding ambient gas into a shell. The material in this shell 
is hot and tries to expand both backwards and forwards, 
giving rise to two shock fronts, one propagating forwards 
into the ambient medium and one backwards into the wind. 
Initially, the shocked wind region can cool efficiently, by in- 
verse Compton scattering for SMBHs (King 2003) and by 
atomic processes for NCs (McLaughlin et al. 2006). As such, 
this region is geometrically thin and the shell is effectively 
driven outwards by a transfer of momentum from the wind 
impacting on its inside. 

The thrust on the shell from the CMC wind is propor- 
tional to the Eddington luminosity of the CMO (King & 
Pounds 2003; McLaughlin et al. 2006): 



dt 



X- 
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where Af^^^Q is the CMO mass and k is the electron scatter- 
ing opacity. For SMBHs, A ~ 1 (King & Pounds); for NCs, 
A ~ 0.05, a value related to the mass fraction of the massive 
stars that contribute to the superwind (McLaughlin et al.). 

As the shell moves outwards, the cooling time of the 
shocked wind material behind the shell eventually becomes 
longer than the dynamical time of the wind. This region 
then cannot cool before more material/energy is injected 
(King 2003; McLaughlin et al. 2006). As such, it expands 
and the shell becomes driven by the thermal pressure in the 
shocked wind region. If the shell can reach a galactocentric 
radius where this switch from momentum- to energy-driving 
occurs, then it may accelerate from that point to escape the 
galaxy (King 2003). 

In this paper, we consider only the momentum- 
conserving phase of the feedback, in the form of a spherical 
supershell moving outwards into a spherical, dark-matter 
dominated protogalaxy, driven entirely by a steady wind 
from a central point mass that may be thought of as either 
an SMBH or an NC. Our aim is primarily to explore the 
effects of relaxing the assumption of isothermal dark mat- 
ter distributions, so we leave to one side all issues around 
any transition to energy-driving, additional feedback from 
bulge-star formation, and evolution of the CMO mass. 

The equation of motion that we consider for the shell is 



^[M,(r).] = A:^-^^ 



[McMo + M^^ir)] , (4) 



where r is the instantaneous radius of the shell; v — dr/dt 
is the velocity of the shell; Mj^^ (r) is the dark matter mass 
inside radius r; and Mg(r) is the ambient gas mass originally 
inside radius r (i.e., the mass that has been swept up into 



the shell when it has radius r). The first term on the right- 
hand side of equation Q is the wind thrust acting on the 
shell, from equation Q . The second and third terms on the 
right-hand side are the gravity of the CMO and the dark 
matter inside the shell (see also King 2005). 

In general, we write Mg(r) = /o h{r) M^j^{r), where /o 
is a fiducial gas fraction (« 0.2) and h{r) is a function that 
describes how the gas traces the dark matter; when h{r) = 1, 
the gas directly traces the dark matter. It is also convenient 
to define characteristic mass and radius scales, Ala- and r^, 
in terms of a characteristic velocity dispersion ao in the dark 
matter halo; 



Ma = /oKcr,t/(A7rG') ~ 4.56 X 10** Mo cr^oo./'o.2A-' 
Ta = GMa/ol ~ 49.25 pcaioo/o.2A-\ (5) 

where CT20o = o-o/200 km s~^ and /0.2 ~ /o/0.2. Referring 
back to equation ([T]), the unit Ma is just the critical CMO 
mass found by King (2005). 

Then, defining M = M/Ma, r = r/va and v = v/gq 
equation Q can be written 

d ^ sM^,,^h{7)M^^{7] 



dr 



h^MlJi\r) 



2h''{7)Ml(T' 

^ ' DM ' 



(6) 



We aim to solve this equation for the velocity fields of 
momentum-driven shells, u ^( r ), rather than r (t) explicitly. 

If the wind thrust is great enough, then equation ([6]) will 
have solutions that allow shells to reach arbitrarily large r 
with non-zero v — the minimum requirement for escape of 
the feedback. If the wind thrust is unable to overcome the 
combined gravity of the CMO and dark matter then the shell 
will stall with « ^ = at some finite radius, and subsequently 
collapse. Equation ((6]) cannot describe such a collapse, since 
that would involve a shell with fixed mass rather than one 
that continually gathers mass [Afg(r) = j^hir) M^^^{r)\ as 
it moves outwards into a galaxy. 

The form of equation ((B]) allows us to select any density 
profile for the host galaxy dark matter and also allows for 
the segregation of gas and dark matter through the function 
h{r). Throughout this paper, we consider only the case that 
h(r) = 1, but we investigate various halo mass distributions. 



3 THE SINGULAR ISOTHERMAL SPHERE 

We look first at the dark matter density profile of a SIS, 
with ^( r ) = 1 so that gas traces the dark matter directly. 
Aspects of this case have been considered previously by sev- 
eral authors (Silk & Rees 1998; Fabian 1999; King 2003, 
2005, 2010a; McLaughlin et al. 2006; Murray et al. 2005). 
King (2005, 2010a) looked at the behaviour of a shell that is 
far from an SMBH, so that the mass of dark matter inside 
the shell dominates over the SMBH gravity. King finds that 
the shell can reach arbitrarily large radii only if the black 
hole has the critical mass given in equation U]). However, 
as we now show, this condition does not actually guarantee 
that a momentum-driven shell will be able to make it to 
large enough radii that the CMO gravity becomes negligible 
compared to the dark matter. 

The density of a SIS is given by 
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singular isothermal sphere 




10-2 0.1 1 10 102 iQZ 0.1 1 10 102 iQZ 0.1 1 10 102 

Figure 1. Velocity fields versus r for momentum-driven shells in a SIS with spatially constant gas fraction and M,-,-^^ = 0.3, 1.01 
and 3. In each case, the solution with C = [r d(0)]^ = is shown by a long-dashed (magenta) line. The physical {v'^ ^0) parts of 
other solutions are shown as solid lines. All solutions with C < are unphysical at small radii, but if M^^^ > 1 they will achieve 
v"^ ^ at large r, corresponding to launches. All solutions with C > decelerate from small radii. If one hits j; = at some point, 
then it generally becomes unphysical at larger radii, and the shell it describes must stall and collapse. This can occur even if M^^^ > 1. 
Formally, solutions with C > and M^j^jq > 1 that stall can have second physical parts with v"^ > at still larger r. These parts of 
such solutions again correspond to launches, though not of the same shells that stall at smaller radii. 



so that 
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In terms of the characteristic mass and radius defined in 
equation ((5]), this means 



Then, with h{r) = 
sheU becomes 

which has solution 



(9) 

1, equation ((Gjl for the motion of the 



2Mr. 



4r 



: 2M^ 



2Mr. 



C_ 



(10) 



(11) 



The constant of integration, C, represents (the square of) 



the shell's momentum, \Mi,(r 



w(r)]^ cx r ^ w 



at r = 0. 



In the limit of very large radius, equation (|lip shows 
that the shell approaches a constant coasting speed: 

~2 



2Mr. 



(r-> 1) 



(12) 



Equations (lllfl and (|12|l are implicit in King (2005) (multi- 
ply his eq. [2] by RR and integrate). Equation (|12p specifi- 
cally is only physical if w ^ > 0. Thus, for the shell to have 
any chance of escaping we must have M^^^ > 1, which is 
exactly the result of King (2005, 2010a). 

In the limit of small radius, the last term of equation 
pif) for V ^ becomes dominant, and the initial momentum 
of the shell (i.e., C) determines the behaviour of the shell. 



If C ^ 0, then « is large and negative at small radii, 
which is unphysical. However, dv^/dr > 0, so it may hap- 
pen that 11 ^ = at some larger radius and increases further 
outwards. The u ^ ^ part of such a solution is physical, and 
the point at which v ^ = Q can be considered as a "launch" 
radius for a (pre-existing) shell initially at rest. 

A launch solution has — Q and dv^ /dr ^ at some 
Haunch- From equation (|ll|l . this requires 



^launch 



1) 



Mr, 



(13) 



Thus, such solutions are only possible for M^^^q > 1, and 
then only starting from radii 



Mr, 



"^launch 



2(Mc 



1 



(14) 



1 ; '^launch 

= 1. 



oo, so launches are not possible 



As McMc 
when _ 

If C > 0, u ^ is large and positive at small radii but 
dv^ /dr < 0, so the shell decelerates but keeps moving out 
into the galaxy, unless and until — Q is reached at some 
finite r. If this happens, then the shell stalls and is not able to 
escape. If = is never realised, then the shell is formally 
able to escape to large radii while purely momentum-driven. 

The stall radius, at which u ^ = 0, is found from equa- 
tion (fTTI) as 



Mr. 



'"stall 



^^^MO - 2C(i\4,,„ - 1) 



(15) 



2(McMo " 1) 

where we have taken the root with the minus sign since this 
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corresponds to the first instance of « ^ = as tlie siiell moves 
outwards. If rstaii is positive and finite, tiie sliell cannot move 
out beyond tliis radius wliile purely momentum driven. 

In tiie case tfiat Af^^Q < 1, f^taa > for any C > 0. We 
can see tliis by noting tliat when M^^q < 1 the discriminant 
in equation (|15|) is always positive and > M^-^^ so that both 
the numerator and denominator are negative, leading to a 
positive rataii- As there are no physical launch solutions and 
the shell always stalls when M^„q < 1, no shell can ever 



escape while purely momentum-driven if M^^^^q 



< 1. 



In the limit that M„ 



1, we find that 



Mr. 



Mr. 



T'stall ~ 



2g(JfcMO-l) 

A/2 
CMO 
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so Tstaii occurs at some positive and finite radius when 
C > and, because riaunch is infinite (equation [14)). when 
-^CMo ~ 1 exactly no shell can escape. 

When rstaii does not exist (formally, when equation [15] 
is complex) , w ^ = is never realised (for solutions with 
C > 0) and the shell is able to reach arbitrarily large radii 
while being purely momentum-driven. This requires 



M„ 



2C{M^ 



1) < 



(17) 



which, for M^j^^ > 1 (as we know this is the only case 
where purely momentum-driven escape is possible), means 
that escape requires 



C> 



Ml 



mo 



1) 



(18) 



If the value of C does not satisfy this constraint, then the 
shell will stall before ever reaching the radii where it could 
coast at the speed given by equation ([12} , even if M^^^^ > 1. 
This is one reason why the critical CMO mass of King (2005, 
2010a) is a necessary but not sujficient condition for the es- 
cape of momentum-driven CMO feedback from an isother- 
mal sphere. 

Figure [T] plots v ^ versus r from equation pi|) for 
M(,„Q = 0.3, 1.01 and 3, with a range of C values in each 
case. The long-dashed (magenta) curve in each panel is the 
solution with C — for that M^^^q . The physical parts of 
solutions with C ^ are shown as solid lines 

The left-hand panel of Figure [1] shows solutions for 
M(,„Q — 0.3. No solution can escape in this case. Those 
with C > all stall at some finite radius (beyond which 
< 0) , while those with C ^ never give physical values 
of > 0. 



Since we know that AI„ 



1 exactly also has no 



escape, the middle panel shows solutions for M^ 



= 1.01. 



In this case, only a few realistic solutions can "escape," and 
those that do tend to a coasting speed of just v ~ 0.14cro 
at large radii (eq. [12]). In order even to reach the radii 
where this applies, shells must have very large velocity at 
small radii (C > 51 from equation [18], which corresponds 
to t; ~ 0.2c at a distance of 1 pc from the CMO if ctq = 
200 km s^^), or else be launched somehow from very large 
radius (riaunch > 50.5 from equation [14], which corresponds 
to r > 2.5 kpc if ao = 200 km s"^). 

Finally, the right-hand panel of Figure [T] illustrates so- 



significant coasting speed of 2ao at large radii, but those 
with C < 9/4 (or v < 15,000 km s"^ at r = 1 pc when 
(To — 200 km s~^) still stall before they are able to make it 
to large radius. Several of the solutions that escape are those 
with a launch radius; these require riaunch > 3/4 (r > 40 pc 
for (TO = 200 km s"^). 

The escape speed from a truncated isothermal sphere 
is Vcsc 2(To at large radius. Our results (equation [12] in 
particular) show that to achieve this escape speed requires 
-^cMo ^ 3. This is another reason why the condition of King 
(2005, 2010a), i.e., simply Af^^^ > 1, is necessary but not 
sufficient for the escape of a purely momentum-driven shell 
from an isothermal sphere. It is also, in essence, the same 
as the objection raised by Silk & Nusser (2010) against ex- 
planations of observed Af^^Q-cr relations as the result of 
outfiows driven by the central objects alone. However, these 
and all other prior results have come from modelling proto- 
galaxies only as SISs. We look now at the effect of allowing 
more realistic descriptions of dark-matter (and ambient gas) 
density profiles. 



4 NON-ISOTHERMAL DARK MATTER 
HALOES 

4.1 General analysis 

Simulated dark matter haloes have density profiles that 
are shallower than that of an isothermal sphere at small 
radii and steeper than isothermal at large radii. Dubinski 
Carlberg (1991) originally fitted haloes with the profile 
of Hernquist (1990), which has p^^^ cc at small radii 
and oc r"* at large radii. The dark-matter profile of 
Navarro, Frenk & White (1996, 1997) also has p^^ oc r"^ 
at small radii, but oc at large radii. Dehnen & 

McLaughlin (2005) develop a family of physically motivated 
halo models that, with pDjj(r) slightly shallower than r~^ 
at small radii and slightly steeper than at large radii, 
match current simulations at least as well as any other fit- 
ting function. 

The circular speed corresponding to all such density 
profiles, K?(r) = GM{r) /r, increases outwards from the cen- 
tre, has a well-defined peak, and then declines towards larger 
radii. This suggests the peak of the circular speed curve as 
a natural point of reference for velocities, radii, and masses 
in realistically non-isothermal haloes. 

We denote the location of the peak in Vc(r) by rpk and 
the value 14 (rpk) = Vc.pk, and we define 



O"0 = 



(19) 



as a characteristic velocity dispersion in order to specify 
unique mass and radius units, M^ and r^, as in equation 
^ above. Then, recaUing that M = M/Ma, r = r/r^, and 
V = v/ao, so that Vc{r) — Mj^-^, (r)/r, we have 



K,pk = 2 
and 

^\MFpk) = Mpk = 2fpk 



(20) 



(21) 



hit ions for M„ 



3. Formally, all solutions now have a 



We now refer all radii to the peak of the circular speed 
curve, defining x = r/rpk; and we introduce a dimensionless 
mass profile, m{x), such that 
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Mj3^, {x) = Mpk m{x) . 
By construction, then, 
M^,,(l)=A7pk ^ m(l) = l . 
Moreover, = Afj^^,, (r)/ r = 2m{x)/x, and thus, 

dlnm^ 



da; 







/ d In m \ 
V din a; ) 



1 



(22) 
(23) 

(24) 



These generic properties of m(x) are important later in our 
analysis (see especially Appendix . 

With these definitions, equation (|6]) for the motion of a 
momentum-driven shell becomes 

[ft^ m'^ u ^ (a;)] = 4M^„q /i(x)m(a;) 



dx 



^ ^CMo h''{x)7n'{x) ^ h'{x)m\x) ^ ^^5) 



pk 



The formal solution of this, when h{x) = 1 for protogalactic 
gas that traces the dark matter directly, is 



m {x)v (a:) = C + 4Afp^jQ / m{u)du 

'0 



du - 



du 



(26) 

'ipk Jo 

where C = m^{0)v^{0) is again a constant of integration 
representing (the square of) the momentum of the shell at 
the origin. 



4.1.1 Velocity fields at small and large radii 

In the limit of small x, we can assume to leading order that 

m{x) — yAx" , (a; < 1, p > 1) (27) 

where p > 1 because we consider only halo density profiles 
that are shallower than isothermal at the centre. Equation 
(|26p then gives 



C_ 



-2p 



(a: < 1, C / 0) 



(28) 



so, as with the SIS, the integration constant, or (the square 
of) the momentum of the shell at r = 0, determines the 
behaviour of the shell. 

If C > 0, then > and dv^/dx < at small radii, 
and the shell decelerates outwards unless and until — 0, 
at which point the shell stalls and then collapses. 

If C < 0, then < at small radii, which is unphys- 
ical; but dv^/dx > 0, so may become positive at some 
non-zero "launch" radius. 

When C = 0, from equation (f^ . 



4M^ 



A P+1 Mpk 2p-l 



(a; < 1, C = 0) 



(29) 



and the behaviour of the shell depends on the specific values 
of A and p, which in turn depend upon the specific choice 
of dark matter density profile. The solution with C — 
corresponds to a shell having zero momentum at a; = and 
is also the value of C that separates initially decelerating 



solutions that either escape or stall (C > 0), from solutions 
that are launched from rest at a non-zero radius (C < 0). 

Since we consider only haloes that are steeper than 
isothermal at large radii, we must have that dlnm/dlnx < 1 
for X > 1. At large radius, the second term from the right- 
hand side of equation (|26p then dominates, so that 



m^{x) 



m{u)du 



{x » 1) (30) 



The velocity field in the limit a: — >■ 00 is therefore completely 
independent of initial conditions (i.e., no C dependence). To 
leading order, v'^ 0{x^~'') with q < 1. Thus, if the sheU 
can make it to large radii at all in a non-isothermal halo 
it must eventually accelerate. This is in contrast to the SIS, 
where a shell at very large radius can only coast at a constant 
speed. It is the steeper-than-isothermal gradient of PD„(r) 
at large radii in realistic dark matter haloes that leads to 
the acceleration. 



4.1.2 Condition for the escape of a particular shell 

Any momentum-driven shell with a velocity field given by 
equation (|26|l . and with C > 0, decelerates as it moves 
outwards from small radii according to equation (|28|) . The 
same is true of shells with C = 0, if the small- a; value of 



V from equation (|29p is positive. Some shells with C < 
and relatively small launch radii can also have w ^ > and 
dv^ /dx < over some range of radius (see below). Mean- 
while, any solution to equation ((25} [with h{x) = 1] accel- 
erates at large radii according to equation (|30|l . Therefore, 
there is a large class of solutions that go through local min- 
ima in u at intermediate radii. We want to know the CMO 
mass required for a particular shell in this class to escape a 
given galaxy. (The only solutions not in this class are some, 
with C < 0, which are launched from large enough radii that 
they only accelerate outwards, and so always escape.) 

If a local minimum in u ^ exists for a particular solution, 
we denote the radius where it occurs by Xmin, and the value 
of the minimum by v ^j^. Putting h{x) = 1 in equation p5|l 



and setting dv /dx = ai x = 



= 4A/^ 



we then obtain 



din x„ 



^^(a^min) 



^ (a^min ) 



(31) 



pk 



If a shell with a given initial momentum (value of C) is 
to escape a dark-matter halo with given m{x), Afpi;, and 
(Jo = Vc,pk/v'2, then we must have w min ^ so the shell does 
not stall (i.e., cross = 0) before it can start accelerating 
outwards. We refer to the case that i; J^j^ = exactly as the 
critical case, and we denote the values of Af^^^ and a;min in 
this case by Afcrit and x^it. Then, from equation (|3ip . 

, -I -1 



Afcrit = 



m (Xcrit) 



1 m(a:crit) 



Also, setting x — x^it 



pk crit 

= 0, and Af„ 



(32) 



it m equa- 



tion (1261) . and using equation (|32p to eliminate Afpk, yields 
Jo°"' [fnix^rit) - m{u)] [m{u)/uf du + C/4 



Af, 



Jo""' [a;crit/'^(3^crit) - ?tVm(ii)] [m{u)/ufdu 
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(33) 



Equating the right-hand sides of equation (|32p and (|33p 
allows us to solve for a^crit , and then Merit , in terms of C and 
the dark-matter halo parameters. The necessary condition 
for the escape of a purely momentum-driven shell with a 
particular value of C is just M^^^ > Merit- 



Equation (|32|l can give a sensible (positive) value for 
Merit only for shell-and-dark matter combinations such that 
Mpk > m(a;crit)/2;crif This is not a problem in general. Mpk 
is the dark matter mass inside the peak of the dark-matter 
circular-speed curve, in units of Mo- — 4.6 x 10* Mq (t|oo 
(equation [S]), and so will be a large number in real galaxies. 
Meanwhile, the function m{x)/x'^ is always equal to 1 at 
X — 1 (equation [23]), so that having A/pk > m(xerit)/2;crit 
at some reasonable value of Xerit is usually assured. 

The density profiles of realistic dark-matter halo models 
are such that d\nm/ d\nx < 2 in the main, the only excep- 
tion being in the very innermost regions of some models (see 
below). Thus, for most values of a^erit, the integral in the de- 
nominator of equation (|33p is positive; while the integral in 
the numerator is always positive. Therefore, this equation 
implies Merit > for any shell with C ^ 0. Launch solutions 
with modest C < can also have Merit > 0, so long as the 
numerator in equation p3p is still positive. If C is too large 
and negative, then formally Merit < 0, which means that 
such solutions do not actually go through minima in « ^ . 
These correspond to shells, launched from large radii, which 
accelerate monotonically outwards to escape regardless of 
the CMO mass. _ 

Below, we will find the necessary Merit for shells that 
have C = Q (i-e., zero momentum at zero radius) in some 
specific dark-matter haloes. We emphasize, however, that 
this is not the only physically meaningful solution. Solutions 
with C > would describe shells that receive an impulse at 
the centre. Solutions with C < could be of interest for 
shells that stall at some radius inside a galaxy during an 
early phase of CMO growth, and are later "re-launched" by 
feedback from the CMO when it is more massive. 



4-1.3 Sufficient condition for the escape of any shell 

Momentum-driven shells with different initial conditions (C 
values) have different values of a;erit and Aferit, given by 
equations H32p and p3p . To compare these values between 
different shell solutions, we differentiate equation (|32p with 
respect to Xerit, for a fixed dark matter mass Mpj,: 



rfMerit 
rfa;erit 



2m^(a::erit) 2;erit 



m(Xerit)/Mpk 



dlnm(a;erit) 

d In Xerit 



1 



1 dm(a;crit) 



2 Afpk ^crit rfXcrit 



.(34) 



By definition, {dlnm/ dlux — 1) = dlnVi^ / dlnx, which is 
positive at a; < 1 and negative for x > 1 (recall equation 
|24)). Hence, dMerit /dxerit > among shells with sufficiently 
small Xerit, and dMerit /rfxerit < among shells with suffi- 
ciently large Xerit. Setting dMerit /dxerit = for a given dark- 
matter mix) and Mpk therefore identifies the momentum- 



driven shell that has the largest critical CMO mass required 
for escape, M"^?-^^ 

^ To find M^^i^, we first solve the equation 
dMerit /rfxerit = for the radius Xe.max at which the 
shell with exactly this critical mass begins to accelerate. 



din m 



din a; 



1 + 



1 dm 



2 Mpk a;e,max dx 

and then use this value in equation (|32p : 



^max 



1 m(Xc,max) 
A'/pk ^c,max 



(35) 



(36) 



The sufficient condition for the escape of any 
momentum-driven shell is simply M^-^^^ J5 Af™j^''. This triv- 
ially includes any launch solutions of the type, mentioned 
above, that do not go through local minima in but only 
ever accelerate outwards. 

We analyze equations (|35p and psp further in Appendix 
m There we show that, in the observationally relevant limit 
of large halo mass Afpk, 



^max 
-iMerit 



1 + 
1 + 



2 A/pk 
1 

Afpk 



d m 
dx^ 



(Mpk » 1) (37) 



to first order in (1/A/pk). That is, in very massive, non- 
isothermal dark matter haloes, the CMO mass that suffices 
to ensure thejBscape of any momentum-driven shell tends 
to the value M^-^^ — > 1; and the radius where the slowest- 
moving shell driven by a CMO with this mass begins to 
accelerate tends to 2;e,max ^ 1, which is the peak of the 
dark-matter circular-speed curve. 



4-1-4 M-a and M-Vc relations 

The dimensional CMO mass that guarantees the escape 
of any momentum-driven shell from a non-isothermal halo 
follows from recalling the definition of the mass unit Mo- 
(equation [5]) and our identification of a characteristic ve- 
locity dispersion in terms of peak circular speed in the 
halo (equation [IS]). For very massive haloes in particular 
(A/pk 3> Mcr), equation (I37p gives approximately 



^max 
J^crit 



M„ 



.fo K 



4 
O"0 



.fo K 



k 

c,pk 



Numerically, 



1.14 X 10* Mq 



/ Ve.pk y 
V200 km s-V 



fo.2 A 



(38) 



(39) 



. 200 km s~ 

in which A ~ 1 describes CMOs that are supermassive black 
holes, while A ~ 0.05 applies for nuclear star clusters (King 
& Pounds 2003; McLaughlin et al. 2006). 

This result reduces to the M-a relation obtained by 
King (2005) for SISs, in which K — v^fo is constant with 
radius. However, there are significant distinctions between 
previous work and our new analysis. 

First, as we have emphasized, M™/^" corresponds in the 
isothermal case to a CMO mass that is necessary but not 
sufficient for momentum-driven feedback to break out of a 
galaxy; while in the non-isothermal case it is a sufficient but 
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not always necessary CMO mass. The CMO masses required 
for the escape of most momentum-driven shells in a given 
dark-matter halo (as obtained from equations [32] and [33] 
above) will be smaller than M'^?^. On these grounds alone, 
theoretical M^f^-ao or M^j^'^-V'cpk relations from our work 
are expected to be something of upper limits to observed 
M-a or M-Vc relations. 

Second, our general treatment shows that the value of 
M^/J" in equation p8[) or p9[l applies only in the limit of 
very large dark-matter halo mass, A/pk 3> Ma-. The exact 
value of the sufficient CMO mass M™/^" for a specific value 
of Mpk must be obtained from equations p5[) and (|36p . 

Third, our results are the first to incorporate explic- 
itly and rigorously the peak value of a dark-matter circular- 
speed curve (and associated velocity dispersion), which is 
a well-defined quantity in any realistic non-isothermal halo. 
This provides a new basis from which to begin addressing 
observational claims of correlations between CMO masses 
and dark-matter halo properties (e.g., Volonteri et al. 2011; 
Ferrarese 2002). 

In the rest of this section, we illustrate the general re- 
sults we have obtained, by looking in detail at their appli- 
cation to three specific dark-matter halo models. 



4.2 Hernquist model haloes 

The first non-isothermal density profile we consider is that 
of Hernquist (1990). This model has been used to fit dark- 
matter haloes from N-body simulations (e.g., Dubinski & 
Carlberg 1991) and also has the advantage that, with it, our 
problem remains analytically tractable. 

The density of a Hernquist sphere is given by 

PomM = ^ (-)"(! + -)" , (40) 
2 7rrg Vro/ V roJ 

where ro is a scale radius and A/tot is the total halo mass. 
In terms of the characteristic mass and radius of equation 
([5|, the mass enclosed inside radius r is 

2 



1.5 



r/rp 
1 -I- r/ro 



(41) 



The circular-speed curve for this model, ~ ^/^^(r )/ r, 
peaks at rpk — ro, so that the mass inside this radius is 



Mu 



^DM('^pk) = Mpk^ 

Defining x = r/rpk = r/ro, we therefore write 



(l + x) 



A/pk m(x) 



(42) 



(43) 



With the above definitions and h{x) = 1 again, equation 
(|25p for the motion of a momentum-driven shell has the 
general solution 



Mr. 



A Mr. 



1 + x 



1 + x 
16 



1 + x 



3 Afpk a; 

In the limit that x 
ll28l) we have 



+ 



21n(l + a;) 



(44) 



5 1 + x ■ 16 
0, m{x) — ^ Ax^ , and from equation 
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Figure 2. Solid lines show, as functions of Mpk, the sufficient crit- 
ical CMO mass, M'^^^, that allows the escape of any momentum- 
driven shell from a Hernquist halo (upper panel); and the radius, 
2;c,maxi at which the slowest-moving shell driven by such a CMO 
begins to accelerate to escape (lower panel). Dashed lines show 
the necessary Af^rit and associated x^i-it for the escape of the 
particular solution with C = 0. We show results only for halo 
masses Afpk > 4, above which x^rit and Xcmax are single-valued 
functions of Mp^. 

« ' ~> , (x « 1, C / 0) (45) 

Id 

or, if C = 0, equation (|29|l instead gives 



^CMo 



3 A'/, 



pk 



(s < 1, C = 0) (46) 



Thus, for halo masses A/pk > 4, all shell solutions with C ^ 
decelerate from large, positive at small radii. 

In the large- j: limit, m{x) — >■ 4 and equation pOf) gives 



{x > 1) (47) 



All solutions tend to the same form at large radii, corre- 
sponding to acceleration outwards that is independent of C, 
as we expect from the general discussion in i\4.1\ 

The solid lines in Figure [5] show, as functions of A/pk, 
the sufficient CMO mass, A/™/^^, that provides for the es- 
cape of any momentum-driven shell from a Hernquist halo 
(upper panel); and the radius, Xcmax at which the slowest- 
moving shell begins to accelerate towards larger radii (lower 
panel). These quantities have been calculated from equa- 
tions (|35|l and ((36}, with m{x) defined in equation ((43}. As 
expected on general grounds (see equation [37]; also Ap- 
pendix |X}, A/^j^" — > 1 and Xcmax — > 1 (denoting the peak 
of the circular-speed curve in the halo) for large A/pk 3> 1. 
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Figure 3. Velocity fields v^(x) for Mq^q = 0.3, 1, 3 in a Hernquist dark-matter halo with spatially constant gas fraction and dimension- 
less A'/pi^ = 4000. This corresponds to a roughly Milky Way-sized halo with rpjj fa 50 kpc, M-p]^ 4.7 X 10"'^"'^ A-Zq, and ctq ~ 140 km s~-^. 
The top axis gives the radius in units of ~ 25 pc ct^^q fo.2 where /0.2 = /o/0.2. The magenta curve represents the solution with 
C = for each value of CMO mass illustrated. As in Figure [T] the physical part(s) of each solution are shown by the solid lines. 



The dashed hnes in Figure [2] show the necessary CMO 
mass, Merit, that allows for the escape from a Hernquist 
halo of shells with C = [m^ ii — >■ as x — > 0] specifically; 
and the radius, icrit, at which this particular shell begins 
to accelerate outwards. In this case. Merit and Xcrit have 
been calculated from equations (|32|) and (|33p . with C set to 
zero and m{x) taken from equation (|43|) . Now, the necessary 
Merit — > 0.93 in the limit Mpk — >■ cxj (versus the sujficient 
A^CTit" — >■ 1), and the acceleration begins at a;erit — ^ 1.46 
(just beyond the corresponding radius for M^„q = M^^). 

Parameters that give a reasonable, model-independent 
summary of the circular-speed curve of the Milky Way dark- 
matter halo are rpk « 50 kpc and T4,pk ~ 200 km s~^ (see, 
e.g., Dehnen et al. 2006; McMiUan 2011). Thus, Mpk = 
rpkK%k/G « 4.7 X 10" Mq; ao = Ve,pk/\/2 « 140 km s'^ 
and Ma « 1.1 x 10* Mq fo.2 A"\ so that Mpk « 4300. 

Figure [3] shows the solutions from equation (|44|) for 
-^OMo = 1 and 3 in a Hernquist halo with Mpk = 4000. 
In each panel the dashed (magenta) curve shows the solu- 
tion with C = 0. As in Figure (TJ the physical parts of each 
solution (i.e., those with ^ Q) are shown as solid lines. 

The left panel shows solutions with M,,^^ = 0.3. Most 
of these represent shells that stall and cannot escape. Unlike 
the SIS however, it is possible to have launch solutions when 
Mp^o < 1. Solutions with C < but very close to zero are 
launched from inside rpk and initially accelerate, then reach 
a maximum velocity and decelerate. When M^^^ = 0.3, 
these solutions all stall at a finite radius. Solutions launched 
from outside rpk either correspond to large and negative C or 
are the formal continuations of solutions that stall at smaller 
radii, go into the unphysical u ^ < regime, but then later 
recover to « ^ > 0. All launch solutions of this type accelerate 



monotonically towards larger radii and therefore escape; but, 
for this value of M^-^^ , they all start from infeasibly large 
launch radii of order a; ~ 10 (i.e., r ~ 10 rpk ~ 500 kpc) or 
more. For large enough positive values of C it is possible for 
a shell to escape without stalling (or being launched from 
a large radius) when Mj-.^^^ = 0.3, although this is again a 
formal result that is not physically plausible. The uppermost 
curve in the left panel of Figure [3] shows one solution that 
evidently only requires Merit < 0.3 to escape this halo; but 
it has C > 10, which, given that rpk ~ 50 kpc and ao « 
140 km s~^, corresponds to a shell velocity of ~ lO'^ c at a 
radius of 1 pc. 

The middle panel of Figure |3] shows solutions for 
-^CMo ~ 1- '-'f solutions shown are able to escape 
the halo. However, from equation (|37p . we know that, with 
Mpk = 4000, the CMO mass sufficient to ensure escape is 
actually M™^^ ^ 1.00025. Thus, there are some shells in a 
(narrow) range of C values very close to that stall rather 
than escape; these are simply not shown here. Several solu- 
tions are shown that have C < but close to 0. These are 
launched from inside the peak of the circular-speed curve 
and, though they come very close to stalling, those shown 
manage eventually to accelerate and escape to large radii. 
Several launch solutions with large and negative C are also 
shown, all starting from radii r > rpk and all accelerating to 
escape. The solution with C = exactly is seen to escape 
(as do all solutions above it with C > 0), which is expected 
since our calculations above (see Figure[2)| gave Merit — 0.93 
for this solution. 

The right-hand panel of Figure |3] shows velocity-field 
solutions for M^^g = 3. This is well above the value of the 
sufficient M^^^^ given by equations (|35p and (|36|) above (or, 
approximately, equation [37]). As a result, and in contrast to 
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the SIS, all shells are able to escape and there are no stalls, 
regardless of the initial shell momentum. 

4.3 NFW model haloes 

We next consider the dark matter density profile of Navarro, 
Frenk & White (1996, 1997; NFW), which has 

-2 



(48) 



where is a scale radius and ps is the density at r^. From 
this. 



ln(l + r/rs) - 



r/r 



1 + r/r 



(49) 



and it follows that the circular-speed curve, = 
GM^^{r)/r, peaks at 7^ = rpk/rs ~ 2.16258. Thus, with 
X = r/rpk we have 



(a; = 1) = = Wnr^.p, ln(l +11)- 

and 

ln(l + Tlx) - TZx/{l + Tlx) _ 



n 



l + TZ 



(50) 



M, 



pk- 



ln(l + 7^)-7^/(l + 7^) 



Mp],m{x) 



(51) 



At small radii, then, the dimensionless mass profile tends to 



m{x) 



-— ln(l + 7^)- 



{x < 1) (52) 



l + 7^J 

implying for the velocity field, from equation 1)28^ for C 7^ 0, 

7^ 



7^4 



ln(l + 11) 



l + 7^ 



or from equation (|29|l for C = 



ln(l + 7^) - 



7^ 



l + 7^ 



(a; < 1, C / 0) (53) 



M^„„ 4 



7^4 



3 Mpk 
(a; < 1, C = 0) (54) 

For Mpk > 5.001, v ^ tends to a positive value in the limit of 
small X for C = 0, and then all shells with C ^ decelerate 
from large, positive velocities at small radii. 

In the limit that a; — > 00, the NFW mass profile diverges 
logarithmically, 

-1 



m(x) 



hi(l +7^) 



7^ 



(l + 7^)_ 

[x > 1) (55) 
and, from equation (|30|) . all shell velocities tend to 



4M, 



'\n{Hx)Y 



ln(l + 7^) - 



7^ 



l + 7^ 



[a; ln(7?.a;) — a;] 
(x> 1) 



(56) 



The solid lines in Figure |4] show, as functions of Mpk, 
the sufficient CMO mass that allows for the escape of any 
momentum-driven shell from an NFW halo (upper panel); 
and the radius, a^cmax, at which the slowest-moving shell 
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Figure 4. Solid lines show, as functions of Mpk, the CMO mass 
A/^^J" , which is sufficient for the escape of any momentum-driven 
shell from an NFW halo (upper panel); and the radius 
at which the slowest-moving shell begins to accelerate to escape 
(lower panel). Dashed lines show the necessary values of M^rit, 
and the associated radii Xcriti for the escape of the particular 
solution with C = 0. Results are shown for Afpk > 5, above 
which Xcrit and Xc,max are single-valued functions of A/pk- 



begins to accelerate (lower panel). These are again calcu- 
lated from equations (|35|l and (|36|l . now with m(a;) given 
by equation (f5T|) . In the limit of large A'/pk, M^jJ:" — >• 1 and 
a^c.max — >■ 1 again, just as found for the Hernquist halo in 
Figure [2] and as expected in general from equation (|37|l and 
Appendix 1X1 

The dashed lines in Figure |4] show the critical CMO 
mass that is necessary for the escape from NFW haloes of 
shells with C = specifically; and the radii, a;crit at which 
these particular shells begin to accelerate for a given Mpk. 
In this case. Merit and Xcrit are calculated from equations 
(|32|l and p3[) . In the limit of large Mpk, we have Merit — >■ 
0.94, again slightly smaller than the CMO mass sufficient 
to ensure the escape of any shell. The acceleration begins 
at Xcrit — >■ 1.50, again somewhat larger than a::c,max in the 
sufficient case. 

Given m{x) in equation (|51|l . equation p6|l for the ve- 
locity fields of momentum-driven shells in NFW haloes must 
be evaluated numerically. Figure [S] shows several of the solu- 
tions for dimensionless CMO masses M^-^^ = 0.3, 1 and 3, 
and with Mpk = 4000 for a Milky Way-sized halo (see 
In each panel of this figure, the dashed (magenta) line shows 
the solution with C = 0. As in Figures [1] and |3l the physical 
parts of solutions with C 7^ are shown by solid lines. 

Figure [5] is qualitatively similar to Figure [3] for shells in 
Hernquist (1990) haloes. The left-hand panel, which plots 
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Figure 5. Velocity fields v'^{x) for M^jj^q = 0.3, 1 and 3 in an NFW halo with spatially constant gas fraction and Afpjj = 4000. Radius 
is shown in units of ~ 25 pc o'j^q/o.2A~^ along the top axis, and in units of Tpi^ 50 kpc along the bottom axis. The dashed, 
magenta curve in each panel represents the solution with C = for that value of A/^^^q. The physical part(s) of all other solutions are 
shown as solid lines. 



solutions for a modest M^^q = 0.3, shows all physically 
plausible shells with C ^ stalling and unable to escape 
the halo. (The one such solution shown that is able to es- 
cape, given this CMO mass, has C > 20, corresponding to 
V ~ 10^ c at r — 1 pc.) Solutions with C < include those 
that are launched from within r < rpk, which first accelerate 
but then decelerate and stall; and those launched from out- 
side r > Tpk, which accelerate monotonically outwards and 
always escape, but which all start from large r > 500 kpc. 

The middle panel shows solutions for M^^^ = 1. All 
of those shown escape, including those with C < but near 
zero, which are launched from r < rpk. There are solutions 
within a narrow range of C values near zero that cannot 
escape. These are not shown, but they exist because, given 
that Mpk — 4000 here, the critical CMO mass required for 
the escape of all possible solutions is Af™]^" ~ 1.00025 > 1, 
according to equation (|37p . The solution with C = is able 
to escape, as the CMO mass necessary to expel it from such 
a massive halo is Afcrit — 0.94 < 1 (see Figure |4]). 

The right-hand panel of Figure [5] confirms again that 
all shells escape easily when M(,„q > M^i^^. 



4.4 Dehnen &: McLaughlin model haloes 

Finally, we consider a dark-matter density profile from the 
family developed by Dehnen & McLaughlin (2005). Their 
models are analytical solutions to the spherical Jeans equa- 
tion, which have "pseudo" phase-space density profiles, 
p{r)/a^{r), that are power laws in radius and closely match 
those found in cosmological Af-body simulations. They also 
allow for radially varying anisotropy in the dark-matter ve- 
locity dispersion; and they fit the spherically averaged den- 



sity profiles of simulated haloes as well as, or better than, 
any other fitting function proposed to date. 

The halo model of Dehnen & McLaughlin that is 
isotropic at its centre has the density distribution 

-6 



(57) 



where ro is a scale radius and Aftot is the total halo mass. 
This gives the enclosed mass profile, 

5 



4/9 



1 + (r/r-o)4/9 



(58) 



The circular-speed curve in this case peaks at rpk/ro 



(11/9) 
and 



9/4 



■ we set X = rjr^^ = (9/ll)^/''(r-/ro). Then, 



[x = \) = Afpk = 



V2oy 



A^DMW = Mpk 



(59) 



A//pkm(a:) . (60) 



When X is small, m(x) — > (20/9)^3;'^''^^, and from equation 
(|28|l the momentum-driven shell velocity field for C 7^ 
tends to 



C(|;)'%-"" 



(2: < 1, C / 0) (61) 



or, from equation (|29p if C = 



Mr. 



-11/9 



(a; < 1, C = 0) (62) 



~2 36 /_9_ 
^ 29 V2O/ 

When X is large, m(x) — >■ (20/11)^ and equation (|30p gives 



( 

la 



(a; > 1) (63) 



12 R. C. McQuillin and D. E. McLaughlin 



1.5 



0.5 




1.5 
1 

0.5 h 



~i — I 1 1 1 1 III — I — I 1 1 1 1 1 II — I — I 1 1 1 1 III — I — I 1 1 1 1 II 



halo: Dehnen k McLaughlin 



I I I M ill — I I M ill — I I I M ill — I I I M il 



_i I I 1 1 1 ill 



10 



I I I 1 1 ill I I 1 r iiiil I I 

102 IQ3 



Figure 6. Solid lines show, as functions of Mpk, the CMO mass, 
M™j^, that is sufficient to ensure the escape of any momentum- 
driven shell from a Dehnen & McLaughlin (2005) halo (upper 
panel); and the radius, a;c,maxi a^t which the slowest-moving shell 
begins to accelerate to escape (lower panel). Dashed lines show 
the necessary Merit) and the associated Xcriti for the^scape of 
shells with C = specifically. Results are shown for Mpjj > 10, 

since then Xcrit and Xc.max are single-valued functions of Mpi^. 



The solid lines in Figure [S] show, as functions of Afpk, 
the CMO mass M'^f^ , which is sufficient for the escape of 
any momentum-driven shell from this halo (upper panel); 
and the radius a;c,max at which the slowest-moving shell 
driven by a CMO with the sufficient mass begins to accel- 
erate outwards (lower panel). These are calculated as usual 
from equations (|35p and (|36p . with m(x) in equation (|60p . 
As for the other haloes we have looked at, and as will always 
be true in general, M™;^" — >• 1 and ic.max — ^ 1 for Mpk » 1. 

The dashed lines in Figure [S] show the necessary CMO 
mass Merit, and the radius Xcrit at which acceleration begins, 
for the escape of shells with C = 0, calculated from equa- 
tions (1311) and (1331)- In the limit of large Mpi^, M^it 0.96 
in this case, and a;crit — > 1.53. 

With m{x) in equation (|60p . the solutions to equation 
H25p with h{x) = 1 must again be obtained numerically. 
Figure 13 shows solutions for several shells in a Dehnen & 
McLaughlin (2005) halo with Afpi< = 4000 (again as in ^42)) . 
for each of the CMO masses M^j^^ — 0.3, 1 and 3. The 



solution with C = in each case is shown by a dashed 
(magenta) line, and the physical part(s) of C 7^ solutions 
are drawn as solid lines. 

Figure 13 is similar in all respects to Figures 13] and IS] 
for the other non-isothermal halo models we have exam- 
ined. The left-hand panel of the figure shows again that 



with M(-,„Q < 1, all physically interesting solutions corre- 
spond to shells that stall. Launch solutions with C < that 
escape must start from impractically large r > 500 kpc. 
Solutions with C > require C > 30 to escape, which im- 
plies unphysical shell speeds at small radii (i.e., v > 10® c at 
1 pc). The middle panel of FigureQconfirms that M^^j^^ = 1 
is almost sufficient for the escape of all momentum-driven 
shells; there are a few solutions with a narrow range of 
C values near C = that cannot escape (because in fact 
A^cTt" ~ 1 + 1/Mpii = 1.00025 here), but which are not 
shown. The right-hand panel finally illustrates again how 
any shell, with any initial conditions, can escape the halo 
when Mp^Q exceeds the sufficient M^^^ given by equations 
H35|) and (|36|l in general. 



5 SUMMARY AND DISCUSSION 

We have analyzed the motion of momentum-conserving su- 
pershells driven into isothermal and non-isothermal proto- 
galaxies by steady (time-independent) winds from central 
massive objects (CMOs: either supermassive black holes or 
nuclear star clusters). Our main goal has been to find the 
critical CMO mass that can drive a supershell to escape a 
galaxy, essentially clearing it of ambient gas and stopping 
further CMO growth. Having such a critical CMO mass as 
a function of a characteristic dark-matter halo velocity dis- 
persion then gives a theoretical M^^^-a relation. 

We assumed that the CMO wind thrust is propor- 
tional to M^-^Q (through the Eddington luminosity: King 

6 Pounds 2003; McLaughlin et al. 2006) to obtain a general 
equation of motion for momentum-driven shells (equation 
in] or 12S]) that allows for any dark- matter halo mass pro- 
file and also for the segregation of gas and dark matter. We 
solved this equation for v'^{r), the (square of the) shell ve- 
locity as a function of radius in the CMO's host galaxy, for 
a number of different dark-matter density profiles, though 
only ever considering the case that gas traces dark matter 
directly. This analysis extends and generalizes others in the 
literature, which have only considered dark-matter haloes 
described as SISs, and which have not presented full solu- 
tions for the velocity fields of momentum-driven supershells. 

Since our main aim was to clarify the effect on theo- 
retical M-a relations of relaxing the simplifying assumption 
that CMOs are embedded in singular isothermal dark matter 
haloes, we retained some other simplifications also adopted 
by previous authors. One of these is the assumption that 
the wind driving the CMO feedback is time-independent — 
in essence, that the CMO mass is constant throughout the 
motion of a momentum-conserving supershell. In reality, of 
course, if the CMO is a black hole emitting at the Edding- 
ton limit, then it is also accreting mass at the Eddington 
rate; thus, a wind thrust proportional to M^^,^ must grow 
on the Salpeter timescale of ~4— 5 x lO'^ yr. If the CMO is 
a nuclear star cluster, then the duration and strength of a 
superwind from it is tied to the star-formation history and 
to the main-sequence lifetime of supernova progenitors. 

The other simplification we made was to consider only 
the momentum-driven phase of supershell evolution, ignor- 
ing any eventual transition to the energy-driven regime. Fur- 
ther work is needed to incorporate gas cooling properly into 
a fuller treatment of time- dependent feedback, which will 
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Figure 7. Velocity fields v ^(x) for CMO masses M^^q = 0.3, 1 and 3 in a Dehnen & McLaughlin (2005) dark-matter halo with spatially 
constant gas fraction and Mp]j = 4000. Radius is in units of ~ 25 pc uf^g fo.2 along the top axis, and in units of rp]j Ri 50 kpc 
(for a Milky Way-sized halo) along the bottom axis. The solution with C = is shown by a dashed (magenta) line in each panel. As in 
Figures [T] [S] and [5] the physical part(s) of all other solutions are shown as solid lines. 



also account for the impact of variable CMO masses and 
wind strengths on M-a relations. 

5.1 The singular isothermal sphere 

Revisiting the case of a galaxy modelled as a SIS, we showed 
in 33] that at large radii a momentum-driven shell tends to 
a constant coasting speed given by equation p2|) : 



(r 00, SIS) (64) 



in which (cf. King 2005) 



fo K 4 



4.56 X 10* Mq atoo fo.2 , (65) 



where ctq is the velocity dispersion of the halo and (J200 = 
cro/(200 km s~^); fo « 0.2 is a fiducial gas mass fraction; 
and the parameter A ~ 1 if the CMO is a supermassive 
black hole (SMBH), 01 X ^ 0.05 if the CMO is a nuclear 
star cluster (NC; McLaughlin et al. 2006). This shows that a 
momentum-conserving shell can reach arbitrarily large radii 
in an isothermal sphere, and potentially escape, only if the 
CMO driving the shell has a mass M^j^^ > Ma (so that 
> 0). Otherwise, any shell must stall at some finite ra- 
dius, and subsequently collapse, until the CMO grows in 
mass and drives a stronger wind (see also King 2005). 



The critical M^^^q value in equation (|65p has previously 
been obtained by methods that did not include solving ex- 
plicitly for v^{r) (see King 2003, 2005, 2010a; Murray et 
al. 2005; McLaughhn et al. 2006). By solving for the full ve- 
locity fields v'^{r) of momentum-driven shells, we have shown 
that, while M^^^ ^ Ma is necessary, it is not sufficient to 
guarantee the escape of momentum-driven CMO winds from 
isothermal spheres. 



First, as discussed in ^ AfcMo '^'^'^ initial momen- 
tum of a shell very near a CMO together determine whether 
the shell can reach large enough radii to achieve the asymp- 
totic coasting speed, v^o', if it cannot, then the value of Vao, 
which is determined by M^^q alone, is immaterial. As an 
example, when M;,^^^ = l.OlMa, a shell will stall at a finite 
radius, and re-collapse, unless its launch from the CMO gives 
it an exceedingly fast velocity of t; > 0.2 c (7200 at a radius 
of r ~ 1 pc cr^oo- 

Second, if a shell is to coast at large radii with the nom- 
inal "escape" velocity from an isothermal sphere — that is, 
with > 2a"o — then our work shows that M^^^q > 3Mo- is 
required. This would mean CMO masses almost an order of 
magnitude higher, at a given ctq, than those provided by the 
observed M-a relations for either SMBHs or NCs. This is, in 
essence, the objection raised by Silk & Nusser (2010) to the 
idea that momentum-driven CMO winds are the sole source 
of M-a. However, the objection — and detailed answers to it, 
whether involving additional feedback from bulge-star for- 
mation triggered by the CMO outflow (Silk & Nusser 2010) 
or a transition to energy-conserving evolution at some large 
shell radius (Power et al. 2011; King et al. 2011) — applies 
only if the host galaxy of a CMO is an isothermal sphere. 



5.2 Non-isothermal haloes 

More realistic descriptions of dark-matter haloes have den- 
sity profiles that are shallower at small radii than the r~^ 
profile of an isothermal sphere, and steeper than r~^ at large 
radii. Therefore, they have circular-speed curves, V^ir) = 
GM{r)/r, with well-defined peaks. We showed that, in any 
such non-isothermal halo, any momentum-driven shell must 
begin to accelerate beyond some large radius and will even- 



14 R. C. McQuillin and D. E. McLaughlin 



tually exceed the halo escape velocity, just so long as the 
CMO wind driving the shell can push it to the radius where 
it starts accelerating. We obtained equations that can be 
solved for the critical CMO mass, Merit, required for the es- 
cape of a shell with a given initial momentum in any halo 
with a peaked Vdr) curve fi )4.1.2l equations [32] and [33]). 
We then showed that there is a largest critical CMO mass, 
M^/J", in any such halo. Once a CMO exceeds this mass, 
any momentum-driven shell can escape the halo ( ii4.1.3|l . 
Our equations psp and p6|) allow the calculation of M'^f^ 
in general and provide a sufficient condition for the escape 
of momentum-driven feedback from non-isothermal haloes. 

In this general analysis, a basic mass unit Mo- is defined 
in terms of the peak circular speed in a halo: 



Mo = 



fo 1^ Ktpk 



A7rG2 



1.14 X 10' Mq fo.2 ( ^^■'"^ A 

V200 kms-V 



(66) 



.200 km s"\ 

In the most relevant case that haloes are much more mas- 
sive than Ma, the sufficient condition for the escape of 
momentum-driven feedback is (equation [37]) 



M" 



= Mo 



1 + 



Mo_ 

Mpk 



+ 



Ml 



(Mpk > Mo) (67) 



where Mpk = ^pkK^pk 



/G is the mass of dark matter inside 
the radius where the halo circular speed peaks. For the Milky 
Way, Afpk « 4000 Mo , so this condition is A/^^o ~ Ma to a 
good approximation in intermediate and massive galaxies. 

In a SIS, 14 = \/2cro is constant and, in effect, A/pk = 
CO, so formally M^/J" and Ma reduce to equation (|65p . Al- 
though important differences remain between the isother- 
mal and non-isothermal cases, this suggests that the most 
appropriate single velocity dispersion to use to characterize 
an entire non-isothermal halo, at least in discussions of M-a 
relations, is simply (tq = Vc,pk/\/2. 

We illustrated the application of our general results by 
solving for the velocity fields of momentum-driven shells in 
three specific models of non-isothermal dark-matter haloes 
(Hernquist 1990— ig^] Navarro et al. 1996, 1997—^^ 
Dehnen & McLaughlin 2005 — ^4.4p . We noted that there are 
two main types of «^(r) solutions, corresponding to shells 
that decelerate from small radii close to the CMO (going on 
to accelerate further out if M^-^^ is large enough, or to stall 
at a finite radius if not), and shells that are launched from 
zero velocity at non-zero radii (and may then either stall or 
escape at larger radii). We also saw that the radius at which 
any particular shell starts to accelerate to escape a halo is 
typically within a factor of order unity times the radius at 
which the dark-matter circular speed peaks (which is some 
tens of kpc in a Milky Way-sized halo). 

Since M^-^^ ^ M™]^" « Mo is a sufficient condition 
for the escape of momentum-driven feedback from non- 
isothermal haloes, it generally exceeds the minimally nec- 
essary condition for the escape of any one particular shell. 
In the specific haloes that we looked at, shells with zero ini- 
tial momentum reach large radii and accelerate to escape 
for any M^-^^ > (0.93-0.96) Mo. Different initial conditions 
may enable escape for still (slightly) lower CMO masses. 



momentum-conserving shells in non-isothermal haloes to 
accelerate at large radii — rather than just to coast as in 
isothermal spheres, at potentially sub-escape speeds even 
if A'/(,j^jQ > Mo — effectively answers the main objection of 
Silk & Nusser (2010) to momentum-driven feedback from 
CMO winds as the direct cause of observed M-a relations. 

Again, these results are for time-independent 
winds from CMOs with fixed masses. We have in- 
tegrated v{r) to find r{t) for the C = Q momentum- 
driven shells in each of the non-isothermal haloes 
calculated in §gl3H131 For A4mo = ^^Tt", these 
shells take ~ 3 — 4 x 10* yrs to move from r = to 
r ~ Tpk, from where they can accelerate to escape 
the galaxy. In the case that the CMO is an SMBH, 
this corresponds to ~ 7—8 Salpeter times. Thus, if a 
critical mass black hole were to launch a momentum- 
driven shell from r = Q, the hole would be a factor 



of 



times more massive by the time the shell 



The fact that M„ 



^ A/" 



allows all purely 



escapes. Our expression for the sufficient M^^^ as a 
function of Vc.pk (eq. [67]) would then presumably 
estimate a lower limit to observed SMBH M—o rela- 
tions. However, this apparent difficulty will be mit- 
igated by two effects. 

First, if SMBHs grew from much smaller seeds, then 
even in the case of purely momentum-driven feedback the 
supershells of swept-up ambient gas will have already been 
driven to large radii by the time the black hole reaches the 
critical mass. The question then becomes, for a given mass- 
accretion history, how near to r = rpk is a supershell at the 
time that the black hole attains our critical mass; and can 
the shell subsequently move out to rpk, and start to acceler- 
ate, within less than another Salpeter time? To answer this 
will require solving a fully time-dependent problem includ- 
ing CMO masses and wind thrusts that (in the SMBH case 
at least) increase monotonically with time. Whatever the fi- 
nal result, it is clear that any upwards "correction" to our 
Merit for steady winds and momentum-conserving shells will 
be substantially less than a factor of ~e^~*. 

Second, the time required for a shell to reach a radius at 
which it can accelerate to escape a galaxy will be less than 
any time we derive, whether for steady or time-dependent 
winds, if the shell transitions from momentum-conserving 
to energy-conserving at some radius (say, rtrans), inside rpk. 
This will happen if the cooling time of the shocked gas in 
the shell exceeds the dynamical time of the wind at r = 
''trans (cf. King 2003; McLaughlin et al. 2006). The issue then 
becomes to find the CMO mass at the time when r — rtrans, 
rather than the mass when r — rpk. 

These considerations emphasize the need to include 
both cooling processes and time-dependent winds in future, 
more sophisticated analyses of CMO feedback-regulated 
galaxy formation. Meanwhile, it is worth noting how tan- 
talizingly close the M-a relations contained in our present 
work already are to the observed scalings. 

5.3 Observational implications 

Although subject to the indicated caveats about 
time-dependent winds and pure momentum-driving, 

our results directly predict a relation between SMBH (or 
NC) masses and the dark-matter haloes of their host galax- 
ies, through the peak circular speed of the haloes (equa- 
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tions [66j and [67], for the large- Mpk limit specifically). 
This provides a basis for understanding relations between 
SMBH mass and dark-matter halo mass or asymptotic cir- 
cular speed, which have been claimed (e.g., Ferrarese 2002; 
Volonteri et al. 2011), though also contested (Ho 2007; Ko- 
rmendy & Bender 2011), on empirical grounds. It is impor- 
tant to recognize the physical content of such a CMO-dark 
matter relation, in a feedback context. It does not suggest 
that dark matter in any way feeds the growth of either black 
holes or nuclear star clusters (cf. Kormendy et al. 2011). 
Rather, it reflects the fact that the gravity of a host galaxy, 
which is dominated by its dark matter halo, is what ulti- 
mately determines whether the feedback from a CMO can 
escape. The more familiar M~a relation has the same fun- 
damental interpretation in this picture. 

Making explicit the connection between a theoretical 
halo 14, pk and an observed stellar a, or even an asymptotic 
circular speed in real galaxies (which will include contribu- 
tions from baryons as well as dark matter), is a nontrivial 
task and beyond the scope of our current discussion. We 
simply recall here that the observed relation between SMBH 
mass and the stellar velocity dispersion averaged over one 
effective radius in a sample of early-type galaxies and bulges 
analyzed by Giiltekin et al. (2009) is 

/ „ \ 4.24±0.41 

Mbh ^ (1.32±0.24) X 10» Mq ——^^ ; (68) 

V 200 km s / 

while the relation inferred by Volonteri et al. (2011) between 
SMBH mass and the asymptotic circular speed in a subset 
of the same systems is 

Mbh ==; (2.45±0.80) X 10' Mo (——^-^—^) .(69) 

V 200 km s / 

If we were to associate our 14, pk and ao = Vc,pk/\/2 in non- 
isothermal haloes directly with observational estimates of 
14, a and (Tcff, then we might conclude that the normaliza- 
tions of the predicted M^^^-Vc^pk and M^^^-a relations 
exceed the observed normalizations by factors of « 3-4. This 
point has previously been made, from comparisons only with 
an isothermal-sphere analysis, by King (2010b). 

However, before too much is made of any nor- 
malization offset, or even the caveats associated with 
steady winds and pure momentum-driving, it is cru- 
cial that the correct relationships be worked out in de- 
tail (within specific dark-matter halo models, and account- 
ing properly for the segregation of dark matter and stars) 
between Vc.pk and 14, a, and between ctq = 14,pk/%/2 and 
the stellar (Teff. It is probably also relevant that we (like 
other authors) have worked with the assumption that the 
gas in protogalaxies directly traces the dark matter. The 
consequences of relaxing this assumption remain unclear, 
although our general equation of motion for momentum- 
driven shells (eq. [6] or eq. [25]) offers a way to investigate 
the question. 

Even with these issues, recognizing the non-isothermal 
structure of real galaxies and dark-matter haloes, and work- 
ing in terms of an M^^jQ-14,pk relation, could provide a 
way to extend and unify discussions and analyses to include 
correlations between CMO masses and host-galaxy proper- 
ties in systems with significant rotational support as well 
as (or even instead of) pressure support. This could be of 
particular interest in connection with nuclear star clusters 



in intermediate-mass ellipticals and bulges, and even in very 
late- type Sc/Sd disks. 
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APPENDIX A: THE MAXIMUM CRITICAL 
CMO MASS 



Equation (|35|l in H4.1.3l is a general expression for the radius, 
Scmax, marking the onset of acceleration of the momentum- 
driven shell that has the maximum critical (necessary) CMO 
mass required to escape a non-isothermal dark-matjter halo 
with a given mass profile m(x) and normalization Mpk: 



din m 



din a; 



1 4- 



1 



1 



dm 



2 Mpk a;c,max dx 



(Al) 



Once this is solved for Xc.max, then equation (|36|) gives the 
value of the maximum critical CMO mass for the halo in 
question: 



.max / 



1 



M, 



pk 



(A2) 



In the limit that i\fpk — >■ oo, the second term on the 
right-hand side of equation (lAip tends to zero, so that 



din m 



din a; 



1 



pk 



(A3) 



But m{x) is defined such that (see equations [23] and [24 
dlnm 



m(l) 



and 



din a; 



x dm 
m dx 



(A4) 



so we conclude that a;c,max — > 1 (the peak of the circular 
speed curve) for large halo masses Mpk — >■ oo. We therefore 



2M 



+ 0{x. 



If (A9) 



pk 



Since the limit Xc,max — >■ 1 corresponds to A/pk — >■ oo, 
terms in (a;c,max — 1)/Mpk are of the same order as terms in 
(a;c,max — 1)^ or terms in 1/Mpj.. With this in mind, solving 
equation (|A9p for X: 



a:c,max = 1 + 7-!-^ J— + 0\-^ 
2K M,k \M,2 



c,max as a function of Afpk gives 

(Afpk > 1) (AlO) 



'pk \^"pk^ 

Finally, putting this into equation (|A2[I yields 



jy^max 



1 



M, 



+ 



pk 



(Mpk > 1) (All) 



As discussed further in Sjl] this is the CMO mass that is 
sufficient to ensure the escape of any momentum-driven shell 
in any non-isothermal halo that has a well-defined peak in its 
circular-speed curve. In general, it is larger than the CMO 
mass that is necessary for the escape of any particular shell. 



for large but finite Afpk (which is the observationally relevant 
situation; see the discussion before Figure [3] in ii4.2|) . which 
also means for values of a;c,max close to 1. 
We define 

j2„ 



d m 

dx"^ 



(A5) 

so expanding m{x) in a Taylor series about x = 1 leads to 

(A6) 
(A7) 
(A8) 



m{x) 

dm 
dx 
dlnm 
din a; 

where we have again used the facts (in equation |A4| ') that 
m = 1 and dm/dx = 1 at a; = 1 always. Equation (lAlf) in 
the limit |a;c.max — II <C 1 is then 



x+- m'i {x - 1)^ + 0{x - ly 
1 + m'i {x-l)+ 0{x - 1)^ 
1 + m" (k - 1) + 0{x - if , 



